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Abstract 

We study the implications of phenomenological processes on a two-singlet extension of the Stan- 
dard Model we introduced in a previous work to describe light cold dark matter. We look into 
the rare decays of T and B mesons, most particularly the invisible channels, and study the decay 
channels of the Higgs particle. Preferred regions of the parameter space are indicated, together 
with others that are excluded. Comments in relation to recent Higgs searches and finds at the 
LHC are made. 
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I. INTRODUCTION 



While still elusive, dark matter is believed to contribute about 23% to the energy budget 
of the Universe jl|. We know it should be massive, stable on cosmic time scales and nonrel- 
ativistic when it decouples from the thermal bath in order to be consistent with structure 
formation. Although its mass and spin are not yet known, masses in the range of 5 — 10 
GeV seem to be favored by the direct detection experiments CoGeNT (2|, DAMA |3| and 
CRESST II |4. 

A number of models have been proposed to try to explain the results of these experiments 
0. As it turns out, having a light dark-matter candidate in supersymmetric theories is quite 
challenging. For instance, in mSUGRA, the constraint from WMAP and the bound on the 
pseudo-scalar Higgs mass from LEP give m x o > 50GeV [6|. Also, in the MSSM, a lightest 
supersymmetric particle with a mass around lOGeV and an elastic scattering cross-section 
off a nuclei as large as 10~ 41 cm 2 is needed in order to fit the CoGeNT data, which in turn 
requires a very large tan/3 and a relatively light CP-odd Higgs. However, such a choice of 
parameters leads to a sizable contribution to the branching ratios of some rare decays, which 
then disfavors the scenario of light neutralinos in the context of the MSSM j3| (see also js]). 

In a recent work (9], we proposed a two-singlet extension of the Standard Model as a 
simple model for light cold dark matter. Both scalar fields were Z 2 -symmetric, with one 
undergoing spontaneous symmetry breaking while the other remaining unbroken to ensure 
stability of the dark-matter candidate. We studied the behavior of the model, in particular 
the effects of the dark-matter relic-density constraint and the restrictions from experimental 
direct detection. We concluded that the model was capable of bearing a light dark-matter 
weakly-interacting massive particle (WIMP) in mass regions where other models may find 
difficulties. We should mention in passing that there are scenarios with unstable light Higgs- 
like particles that have been previously studied in certain extensions of the Standard Model, 
see for example 10]. There is also the possibility of having a light pseudo-scalar in the 
NMSSM, see for example |U|. 

The present work studies the effects and restrictions on the two-singlet model coming 
from particle phenomenology. A limited selection of low-energy processes has to be made, 
and we choose to look into the rare decays of T and B mesons. We limit ourselves to small 
dark-matter masses, in the range 0.1 — 10 GeV. We also study the implications of the model 
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on the decay channels of the Higgs particle and make quick comments in relation to recent 
finds at the LHC. 

The theory starts effectively with eight parameters 9|. The spontaneous breaking of 
the electroweak and Z 2 symmetries introduces the two vacuum expectation values v and v i 
respectively. The value of v is fixed experimentally to be 246GeV and we take v\ = lOOGeV. 
Four of the parameters are the three physical masses tuq (dark- matter singlet Sq ), ni\ (the 
second singlet Si) and uih (Higgs h), plus the mixing angle 9 between h and Sj. We let m,i 



vary in the interval 0.1 — 10 GeV and fix the Higgs mass to m^ = 125GeV [12l.ll3j]. except in 
the part about the Higgs decays where we let vary in the interval 100 — 200 GeV 1 . For 
the purpose of our discussions, it is sufficient to let 9 vary in the interval 1° — 40°. The last 
parameters are the two physical mutual coupling constants Aq 4 ^ (dark matter - Higgs) and 
(dark matter - Si particle). In fact, is not free as it is the smallest real and positive 
solution to the dark- matter relic density constraint 9], which is implemented systematically 
throughout this work. Thus we are left with four parameters, namely, tuq , mi, 9 and Aq 4 \ 
To ensure applicability of perturbation theory, the requirement rj^f < 1 is also imposed 
throughout, as well as a choice of rather small values for Aq 4 \ 



II. UPSILON DECAYS 

We start by looking at the constraints on the parameter space of the model coming from 
the decay of the meson T in the state nS (n = 1, 3) into one photon 7 and one particle Si. 
For mi < 8GeV, the branching ratio for this process is given by the relation: 

Br (T nS -> 7 + Si) = G F<^° Xn U _ ^lf( Xn )) BrW 6 (m TnS - mi) . (2.1) 

In this expression, x n = (l - m\/m\ ns ) with m Tl(3)s = 9.46(10.355)GeV the mass of Yi(3)5, 
the branching ratio Bi ip) = Br (T 1(3)s = 2.48(2.18) x 10" 2 Q, a is the QED 

coupling constant, a s = 0.184 the QCD coupling constant at the scale m TnS , the quantity Gp 



is the Fermi coupling constant and rrib the b quark mass 



the effect of QCD radiative corrections given in 17 ]. 



161 ] . The function f(x) incorporates 



1 The exclusion mass range reported by the CMS and ATLAS Collaborations applies to the SM Higgs and 
can be weakened or evaded in models where the Higgs production and/or decay channels are suppressed 



14j . We will comment on this possibility within our model in the last section. 
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However, a rough estimate of the lifetime of S\ indicates that this latter is likely to 
decay inside a typical particle detector, which means we ought to take into account its most 
dominant decay products. We first have a process by which Si decays into a pair of pions, 
with a decay rate given by: 



T {Si -»■ mi) 



GfTXIi 



sin 2 d 



m| 
27 



1 + 



lime 
2mJ 



x 



1 - 
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(2.2) 



Here, m^ is the pion mass and mx the kaon mass. Also, chiral perturbation theory is used 
below the kaon pair production threshold, and the MIT bag model above, with the dressed 
u and d quark masses M u = Md = 0.05GeV. Note that this rate includes all pions, charged 
and neutral. Above the 2rriK threshold, there is the production of both a pair of kaons and 
rj particles. The decay rate for K production is: 



T (Si -»■ KK) 



9 3G F M 2 mi . 2/1 / 4m 2 K 
■ sin 01 7T- 



G (mi - 2m^) . (2.3) 
I9J. For 7] production, replace 



13 4y/2ir V m . 

In the above rate, M s = 0.45GeV is the s quark bag-mass [1 
m K by m v and ^ by ^. 

The particle Si decays also into c and b quarks (mainly c). Including the radiative QCD 
corrections, the corresponding decay rates are given by: 



T(Si -> qq) 



?>GfVii mi 



sm 
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(l + 5-67^) 6 (mi - 2m q ) . 



(2.4) 



4V2tt 

The dressed quark mass m q = m q {mi) and the running strong coupling constant a s = 
Q s (mi) are defined at the energy scale mi {20]]. There is also a decay into a pair of gluons, 
with the rate: 
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(2.5) 



Here, a' s = 0.47 is the QCD coupling constant at the bag-model scale. 

We then have the decay of Si into leptons, the corresponding rate given by: 

/ 4m 2 



r (Si -». e + e~) 
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(2.6) 



where rri£ is the lepton mass. Finally, S\ can decay into a pair of dark matter particles, with 
a decay rate: 



r (5*1 — » SqSq) 



(3) 



Ami 



6 (rrii — 2m ) . 



(2.7) 



327rmi 

The coupling constant 77^ is given in [9|. The branching ratio for Y„s decaying via Si into 
a photon plus X, where X represents any kinematically allowed final state, will be: 



Br (T nS ^ 1 + X) = Bi (T nS ->• 7 + Si) x Br X) 



(2.8) 



In particular, X = So So corresponds to a decay into invisible particles. 
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FIG. 1: Typical branching ratios of T\s decaying into r's, charged pions and charged kaons as 
functions of mi . The corresponding experimental upper bounds are shown. 



The best available experimental upper bounds on lS'-state branching ratios are: (i 



Br (Tig -»■ 7 + rr) < 5 x 1(T 5 for 3.5GeV < m x < 9.2GeV |2l|; (ii) Br (T 15 7 + tt+tt-) < 
6.3 x lO" 5 for lGeV < m 1 22]; (hi) Br (T iS 7 + X+X-) < 1.14 x 10~ 5 for 2GeV < mi < 
3GeV [23|]. Figure [U displays the corresponding branching ratios of T15 decays via S% as 
functions of m 1; together with these upper bounds. Also, the best available experimental 
upper bounds on T35 branching ratios are: (i) Br (T 3S — > 7 + ////) < 3 x for lGeV < 
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mi < lOGeV; (ii) Br (T 3S -> 7 + Invisible) < 3 x 1CT 6 for lGeV < m 1 < 7.8GeV [24 1. 
Typical corresponding branching ratios are shown in figure |2j 



A (4) = 0.02, m = 1.5GeV, = 20° 




1 2 3 4 5 6 7 8 

mi (GeV) 



FIG. 2: Typical branching ratios of T35 decaying into muons and dark matter as functions of mi. 
The corresponding experimental upper bounds are shown. 

A systematic scan of the parameter space indicates that the main effect of the Higgs- 
dark-matter coupling constant Aq 4 '' and the dark-matter mass m is to exclude, via the relic 
density and perturbativity constraints, regions of applicability of the model. This is shows in 
figures [TJ and [2] where the region mi < 1.4GeV is excluded. Otherwise, these two parameters 
have little effect on the shapes of the branching ratios themselves. The onset of the SqSq 
channel for mi > 2mo abates sharply the other channels and this one becomes dominant by 
far. The effect of the mixing angle 9 is to enhance all branching ratios as it increases, due to 
the factor sin 2 9. The dark matter decay channel reaches the invisible upper bound already 
for 9 ~ 15°, for fairly small mo, say 0.5GeV. The other channels find it hard to get to their 
respective experimental upper bounds, even for large values of 9. 
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III. B MESON DECAYS 



Next we look at the flavor changing process in which the meson B + decays into a K + 
plus invisibles. The corresponding Standard-Model mode is a decay into K + and a pair of 
neutrinos, with a branching ratio Br M (B + — > K + + vv) ~ 4.7x 10~ 6 [25]. The experimental 
upper bound is Br Exp (B + — > K + + Inv) ~ 14 x 10~ 6 [261] . As in the T decays, the most 
prominent B invisible decay in this model is into SqSq via Si . The process B + — > K + + Si 
has a the branching ratio: 

Bl{B ^ K +Sl) ~ 1024^ K-m s ) 2 {VM /o(mi) 



x y {m\ — m\) {m 2 _ — m\) sin 9 (m_ — mi) . (3.1) 

Here m± = tub ± mx where is the B + mass, tb its lifetime, andVj^ and Vt s are flavor 
changing CKM coefficients. The function /o (s) is given by the relation: 



fo (s) = 0.33 exp 



0.63s 0.095s 2 0.591s 3 



+ 

m B m'g m# 



(3.2) 



The different Si decay modes are given in (12. 2p - ( 12.71) above. The branching ratio of B + 
decaying into K + + SoSo via the production and propagation of an intermediary Si will be: 

Br (Si) ( 5 + _^ K + + SoS ^ = Br ( B + _^ K + + s ^ x Br ( Si _^ SqS ^ _ (3 3) 

Figure |3] displays a typical behavior of Br (Sl) (B+ K + + S S ) SIS db function of vri\. 
The branching ratio is well above the experimental upper bound, and 9 as small as 1° will 
not help with this, no matter what the values for Aq 4 '* and mo are. So, for mi < 4.8GeV, this 
process excludes the two-singlet model for m < mj/2. For m x > 4.8GeV or m > mi/2, 
the decay does not occur, so no constraints on the model from this process. 

Another process involving B mesons is the decay of B s into predominately a pair of 
muons. The Standard Model branching ratio for this process is Br SM (B s — > = (3.2 ± 

0.2) x 10 -9 [27], and the experimental upper bound is Br Exp (B s ^^~) < 1.08 x lO" 8 
281 ] . In the present model, two additional decay diagrams occur, both via intermediary Si, 
yielding together the branching ratio: 

x 9T B Gif 2 R m% „ l2 (l-4m 2 /m|) 3/2 , , s 
Br<*> (B s ->^V = 9n ^! s 5 Bs ™lmt\V tb V t :\ 2 -±- gJ sin 4 g. (3.4) 

20487T " ( m 2 _ m 2) + m 2p2 
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A (4) = 0.01, m = 0.2GeV, 6» = 10° 
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FIG. 3: Typical branching ratio of i? + decaying into dark matter via S\ as a function of m\. The 
SM and experimental bounds are shown. 

In this relation, t# s is the B s life-time, mg s = 5.37GeV its mass, and fg a a form factor that 
we take equal to 0.21GeV. The quantity Ti is the total width of the particle Si [91]. 

A typical behavior of Bt^ Si \B s — > ^ //~) as a function of mi is shown in figure HI The 
peak is at m^,. All three parameters A^ , mo and 9 combine in the relic density constraint 
to exclude regions of applicability of the model. For example, for the values of figure HI the 
region mi < 2.2 GeV is excluded. However, a systematic scan of the parameter space shows 
that outside the relic density constraint, Aq has no significant direct effect on the shape of 
3t^ Si \B s — > As m increases, it sharpens the peak of the curve while pushing it 

up. This works until about 2.7GeV, beyond which mo ceases to have any significant direct 
effect. Increasing 9 enhances the values of the branching ratio without affecting the width. 

Also, for all the range of mi, all of Br (5l) + Br SM stays below Br Exp as long as 9 < 10°. 
As 9 increases beyond this value, the peak region pushes up increasingly above Br Exp and 
thus gets excluded. Hoping for a clear signal if any, figure [5] displays a density plot showing 
Br SM +Br^(5 s — > ) in the plane (mi, 9), squeezed between Br M +5cr from below and 
Br Exp from above. The behavior we see in this figure is generic across the ranges of m and 



S 



A (4) = 0.05, m = 3GeV, = 15° 




mi(GeV) 

FIG. 4: Typical behavior of Br^ Sl \B s ) as a function of mj, together with the SM and 

experimental bounds. 

Aq 4 '*: the V-shape structure in gray developing from mi = rriB s is the allowed region. The 
white region in the middle is due to Br Exp , and the white region outside to Br SM + 5a. It can 
happen that some of the gray V is eaten up by the relic-density constraint and perturbativity 
requirement for larger values of Aq 4 ''. 

From this process, there is probably one element to retain if we want the model to 
contribute a distinct signal to B s — > fi + fi~ for the range of mo chosen, and that is to restrict 
4GeV < m 1 < 6.5GeV. No additional constraint on m is necessary while keeping < 0.1 
to avoid systematic exclusion from direct detection is safe. 

IV. HIGGS DECAYS 

We finally examine the implications of the model on the Higgs different decay modes. In 
this part of the work, we allow the Higgs mass to vary in the interval lOOGeV — 200GeV. 
First, h can decay into a pair of leptons £, predominantly r's. The corresponding decay 
rate F (h — > £ + £~) is given by the relation f)2.6p where we replace mi by m^ . It can also 
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V 4) = 0.03, m = 5GeV 
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FIG. 5: A density plot of B SM + Br^ 1 ^ squeezed between B SM + 5a from below and B Exp from 
above. 

decay into a pair of quarks q, mainly into fe's and, to a lesser degree, into c's. Here too the 
decay rate T(h — > qq) is given in (12. 4p with the replacement instead of mi. Then the 
Higgs can decay into a pair of gluons. Including the next-to-next-to-leading QCD radiative 
corrections, the corresponding decay rate can write like this: 



where the sum is over all quark flavors q. A systematic study of the double integral above 
shows that, with in the range lOOGeV - 200GeV, the t quark dominates in the sum over 
q, with non- negligible contributions from the c and b quarks. 

For rrih smaller than the W or Z pair-production threshold, the Higgs can decay into a 
pair of one real and one virtual gauge bosons, with rates given by: 
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r (h vv*) 



3G 2 F mym h 



cos 2 9 A v R I — j ] [(m h - m v ) (2m v - m h )\ . 




(4.2) 



167T 3 
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In this expression, my is the mass of the gauge boson V, the factor Ay = 1 for W and 
(X — y sin 2 TO + y s i n4 f° r % with W the Weinberg angle, and we have: 

3(l-8x + 20x 2 ) /3a; -1 

R[x) = . arccos 



1 — x 3 

(2 - 13a; + 47x 2 ) - - (l - Qx + Ax 2 ) logx. (4.3) 



2x v y 2 

For a heavier Higgs particle, the decay rates into a V pair is given by: 



Gpiriy cos 9 f Am 



1 Q™2 \2 



1 + 



(mi — 2m 



v) 



6 (m h - 2my) , (4.4) 



r{h -> ^ " f, 7""" " ^y ( i - — y 

y2^m h V m 

with 5y = 1 for and | for Z. 

While all these decay modes already exist within the Standard Model, the two-singlet 
extension introduces two additional (invisible) modes, namely a decay into a pair of So's 
and a pair of Si's. The corresponding decay rates are: 



Af A _ 4m: 
327rmh \ m/ h 



r(h^S l S z ) = ^±-(l-—f) e(m h -2m t ), (4.5) 



where A, = A^L f° r So(i) are coupling constants given in [9(, functions of the param- 
eters of the theory. The total decay rate T (h) of the Higgs particle is the sum of 
these partial rates. The branching ratio corresponding to a particular decay will be 
Br (h -> X) = T (h ->■ X) /r (/i). 

Typical behaviors of the most prominent branching ratios are displayed in figure [6j A 
systematic study shows that for all ranges of the parameters, the Higgs decays dominantly 
into invisible. The production of fermions and gluons is comparatively marginal, whereas 
that of W and Z pairs takes relative importance towards and above the corresponding 
thresholds, and more significantly at larger values of the mixing angle 9. 

However, the decay distribution between So and Si is not even. The most dramatic effect 
comes from the coupling constant A 4 \ When it is very small, the dominant production is 
that of a pair of Si. This is exhibited in figure El for which \ { 4) = 0.01. As it increases, there 
is a gradual shift towards a more dominating dark-matter pair production, a shift competed 
against by an increase in 9. Figure [7] displays the branching ratios for Ag 4 ' ) =0.1 and figure 
[8] for the larger value A 4 ^ = 0.7. In general, increasing 9 smoothens the crossings of the 
WW and ZZ thresholds, and lowers the production of everything except that of a pair of 
Si, which is instead increased. 
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FIG. 6: Branching ratios for Higgs decays. Very small dark-matter Higgs coupling. 

Like in the Standard Model, the production of a pair of b quarks dominates over the 
production of the other fermions, and all fermions are not favored by increasing A 4) . Chan gcs 
in mo and mi have very little direct effects on all the branching ratios except that of SqSq 
production where, at small 9, increasing mi (mo) increases (decreases) the branching ratio, 
with reversed effects at larger 9. Note though that these masses have indirect impact through 
the relic density constraint by excluding certain regions {9]. 

V. CONCLUDING REMARKS 

In this work, we have explored some phenomenological aspects of a two-singlet extension 
of the Standard Model we proposed as a simple model for light cold dark matter. We have 
looked into the rare decays of T and B mesons and studied the implications of the model 
on the decay channels of the Higgs particle. 

For both T and B decays, apart from combining with the other two parameters in the 
relic-density and perturbativity constraints to exclude regions of applicability of the model, 
the Higgs-dark-matter coupling constant Aq 4 "* and the dark-matter mass tuq have little effect 
on the shapes of the branching ratios. Also, the effect of increasing the h — Si mixing 
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A (4) = 0.1, 6 = 10°, m = 2GeV, m\ = 2.5GeV 
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FIG. 7: Branching ratios for Higgs decays. Small dark-matter Higgs coupling. 

angle 6 is to enhance all branching ratios. For T decays, the dark-matter channel, when 
kinematically allowed (mi > 2mo), dominates over the other decay modes. It reaches the 
experimental invisible upper bound for already fairly small values of 9 and m®. From B + 
decays, we learn that our model is excluded for m\ < 4.8GeV (= mg-mjf) and m < m\/2. 
From B s decay into muons, we learn that for the model to contribute a distinct signal to 
this process, it is best to restrict 4GeV < mi < 6.5GeV with no additional constraint on 
mo- Also, in general, keeping Ag 4 '' < 0.1 to avoid systematic exclusion from direct detection 
for all these processes is safe. 

Before closing this section, we comment on the effect that dark matter in our model has 
on Higgs searches. Since m h ^> 2m , the process h —> S So is kinematically allowed and, for 
a large range of the parameter space, the ratio 

(6) _ Br (h -> SqSq) 

can be larger than one for < 120GeV as can be seen in figure In this situation, the 
LEP bound on the Higgs mass can be weaker. Also, in our model, the Higgs production 
at LEP via Higgstrahlung can be smaller than the one in the Standard Model, and so the 
Higgs can be as light as lOOGeV. Such a light Higgs would be in good agreement with the 
electroweak precision tests . 
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A (4) = 0.7, 6 = 10°, m Q = 2GeV, mi = 2.5GeV 
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FIG. 8: Branching ratios for Higgs decays. Larger dark-matter Higgs coupling. 
As to the Higgs searches at the LHC, the ATLAS and CMS collaborations have reported 



recently the exclusion of a Higgs mass in the interval 145 - 460 GeV 12, ll3] , which seems 
to suggest that we should have limited our analysis of the Higgs branching ratios to rrih < 
145GeV. However, it is important to note that these experimental constraints apply to the 
SM Higgs and can not therefore be used as such if the Higgs interactions are modified. In 
our model, the mixing of h with Si will result in a reduction of the statistical significance 
of the Higgs discovery at the LHC. Indeed, the relevant quantity that allows one to use the 
experimental limits on Higgs searches to derive constraints on the parameters of the model 
is the ratio: 

v = a (99 -> h) Br (h -S> X SM ) cos 4 6 , , 

XsM ~ <r(SM) ( gg _> h) Br SM (h -> X SM ) ~ cos 2 9 + r (h -> X inv ) /r| M " 1 ' } 

In this expression, A SM corresponds to all the Standard Model particles, X inv = S S and 

S\Si, a is a cross-section, Br SM (h — > X) the branching fraction of the SM Higgs decaying 

into any kinematically allowed mode X, and the total Higgs decay rate in the Standard 

Model. To open up the region m h > 140GeV requires the ratio 7^x SM t° be smaller than 0.25 



131 ]. a constraint easily fulfilled in our model. By comparison, the minimal extensions 



of the Standard Model with just one singlet sea 
symmetry, are highly constrained in this regard 
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ar or a Major ana fermion, even under Z2 



29|]. 



Beyond the present investigation, many other aspects of the model, partly in relation 
with the current LHC searches, are just waiting to be explored. 
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